


SPECIFICATION 



ME I'HOD FOR MEASURING DIGGING POSITION 




BACKGROUND OF THE INVENTION 
[Technical Field of tlie bivention] 

The present invention relates to a method for determining the digging 
position in a non-open-cut metliod of excavation and, more particularly, to a 
5 method which ensures iaccuracy in detemiining positions by lessening tlie 
influence of a noise magnetic field of frequency components close to tliat of the 
signal magnetic field to be measured. 

[Prior Art] 

A horizontal drilling method, which is one of the non-opeh-cut rnethods of 
10 this kind, uses a small-diameter pipe of 100 mm or less across tor horizontally 
digging in tlie ground, and accordingly, such a precision position determining 
apparatus as used in an ordinary small-diameter driving metliod of excavation 
cannot be placed near a drill. To solve tliis problem, it is customairy in the art to 
adopt a method in which an AC magnetic field is generated by a coil mounted in 
15 the drill head and detected by an above-ground magnetic sensor like a coil to 
determine the current digging position. 

Tliis metliod is simple and easy, but since tiie magnetic field by the coil is 
a dipole magnetic field, it rapidly attenuates witli distance from tlie coil. Hence, 
this method has a defect of inability to achieve high-reliability determination of the 
20 digging position when a power line or similar magnetic noise source is present in 
the vicinity of the place where to perform the position determination. 

SUMJMAR OF THE INVENTION 
All object of tlie present invention is to provide a position determination 
metliod tliat pemiits high-reliability deteraiination of the digging position by 
25 detecting on the ground the AC signal magnetic field provided from a coil housed 
in the drill head even if a noise magnetic field is present which affects the position 
deteraiination. 

To attain tlie above object, a digging position determining method 
according to tlie present invention for non-open-cut excavation, which senses an 




AC magnetic field provided from a magnetic field source by a magnetic sensor 
provided on tlie ground and calculates the position of the magnetic field source 
from tlie magnitude and direction of die sensed magnetic field, said metliod having 
a construction characterized in that 
5 In case where in addition to a signal magnetic field generated by said 

magnetic field source, there exists a noise magnetic field generated by a nearby 
current, 

at lease one of die position of said magnetic field source, the tilt angle of 
said magnetic field source to the vertical direction and the azimuth angle of its 

10 axial direction of said magnetic field source in a horizontal plane is calculated, 
from a projective component of die magnetic field sensed by said magnetic sensor 
and projected on a plane or straight line ortliogoiial to a vector-valued direction of 
said noise magnetic field. 

That is, ill 1999 year's investigation and research relating to usefiil 

15 utilization techniques of energy resources: entitled "Research for low-loss 
energizing techniques in establishment of advanced telecommunication 
network",(executed by Composite Development System for New Energy 
hidustrial Technique), it is described tliat extemal noise magnetic fields, which 
affects the position determination in the non-open-cut method of excavation, is 

20 mostly generated by a current of some kind, hi this case, the magnitude of the 
noise magnetic field varies irregularly with time, but its vector-valued direction is 
constant at each field sensing position. 

The present invention attains its object through the adoption of tlie 
following steps (A) and (B). 

25 (A) The direction of the noise magnetic field is detected, and a sensed 

magnetic field in wliich die noise magnetic field and a signal magnetic field are 
mixed is projected on a plane or straight line ortliogonal to the direction of die 
noise magnetic field to obtain a projective component. 

Since tiie projective component is theoretically free from a component 




derived from tlie noise magnetic field, at lease one of the position, azimuth angle 
and tilt angle of the magnetic field source is calculated so tliat tlie magnitude of the 
projective component (in Uie case of projection on tlie straiglit line) or it magnitude 
and direction (in the case of projection on the plane) is substantially equal to a 

5 theoretically calculated value of a corresponding quantity of a magnetic field 
generated from a magnetic field source, or ,has a minimum difference between tlie 
former and tlie latter The sensed magnetic field is obtained by sensing magnetic 
fields at different positions whose number is determined by how many ones of tlie 
position, azimuth angle and tilt angle of the magnetic field source are unknowns 

1 0 and how such unknowns are ceilculated. 

(B) To obtain the noise magnetic field, 

a) Wlien tlie number oF noise magnetic fields is virtually one: 

• The noise magnetic field is sensed to obtain its direction essentially in the 
absence of a signal magnetic field. 
15 . When the noise magnetic field has frequency components at frequencies 

different from diat of the signal magnetic field, the frequency components are 
measured to obtain the direction of the noise magnetic field. 

b) Wlien the number of magnetic fields is virtually two: 

The frequency components of a first one of the two noise magnetic fields, 
20 which are widely spaced from the frequency compoiients of the second noise 
magnetic field and tlie signal magnetic field, are measured to obtain a vector 
direction of the first noise magnetic field, and the frequency components of tlie 
second noise magnetic field, which are widely spaced from the frequency 
components of the first noise magnetic field and the signal magnetic field, are 
25 measured to obtain a vector direction of the second noise magnetic field. 

BRIEF DESCRIPTION OF THE DRAWINGS 
The features of tlie present invention will be clearly understood from tlie 
following description taken in conjunction with the accompanying drawings, in 
which: 




Fig. I is a perspective view explanatory of the placement of a magnetic 
sensor in the present invention metliod; 

Fig, 2 is a vector diagram explanatory of the principle of measurement 
according to the present invention method in a case of one noise magnetic field 
5 source; 

Fig. 3 is a vector diagram explanatory of the principle of measurement 
according to the present invention metliod in a case of two noise magnetic field 
sources; 

Fig. 4 is a flowchart illustrating a procedure for obtaining a projective 
10 component according to the present invention method in a case of one noise 
magnetic field source; 

Fig. 5 is a flowchai t illustrating a procedure for obtaining a projective 
component according to die present invention method in a case of two noise 
magnetic field sources; 
1 5 Fig. 6 is a flowchart illustrating a procedure for obtaining a signal magnetic 

field in a case of one noise magnetic field; 

Fig. 7 is a flowchart illustrating a procedure for obtaining a signal magnetic 
field in a case of two noise magnetic fields; 

Fig. 8 is a flowchart showing a measurement procedure in the present 
20 invention metliod; 

Fig. 9 is a flowchart showing a procedure for calculating the position of the 
signal magnetic field according to the present invention metliod in a case of one 
noise magnetic field when tlie number of unknowns and Ihe number of equations 
are the same as to each other; 
25 Fig. 10 is a flowchart showing a procedure for calculating the position of 

tlie sigiiiai magnetic field source according to the present invention metliod in a 
case of one noise magnetic field source when the number of equations is larger 
than the number of unknowns; 

Fig. 11 is a flowchart showing a procedure for calculating the position of 




the signal magnetic field according to the present invention method in a case of 
two noise magnetic field sources when the number of unknowns and tlie number 
of equations are the same as to each otlier; 

Fig. 12 is a flowchart showing tlie procedure for calculating the position of 
5 tlie signal magnetic field source according to tlie present invention method in a 
case of two noise magnetic field sources when tlie number of equations is larger 
than die number of unknowns; 

Fig. 13 is a flowchart showing an example of a procedure for obtaining the 
signal magnetic field tlirough tlie use of the magnitude of the signal magnetic field 
1 0 vector; 

Fig. 14 is a flowchart showing another example of a procedure for 
obtaining tlie signal magnetic field tlirough the use of the magnitude of the signal 
magnetic field vector; 

Fig. 15 is a perspective view explanatory of the placement of a magnetic 
15 sensor for detemiining the direction of a noise magnetic field according to die 
present invention method when a digging head, which is a signal magnetic field 
source, is distant fi^oni the position of measurement; 

Fig. 16 is a flowchart showing a procedure for calculating tiiie direction of 
the noise magnetic field according to the present invention method in a case where 
20 tlie digging head is distance fi-om the position of measurement and the signal 
magnetic field source stops to generate die magnetic field; 

Fig. 1 7 is a flowchart showing a procedure for calculating the direction of 
the noise magnetic field according to die present invention method when the signal 
magnetic field and the noise magnetic field are mixed to each other; 
25 Fig. 18 is a signal frequency spectrum diagram explanatory of how the 

operation of selecting a frequency for maximizing a frequency spectrum is used in 
the process for calculating the direction of die noise magnetic field in die flowchart 
shown in Fig. 13; 

Fig. 19 is a flowchart explanatory of a first method by which the flow of 



candidate vector calculating process is used in tlie process for calculating tlie 
direction of tlie noise magnetic field in tlie flowchart shown in Fig. 17; 

Fig. 20 is a flowchart explanatory of a second method by which tlie flow of 
candidate vector calculating process is used in the process for calculating the 
5 direction of the noise magnetic field in tlie flowchart shown in Fig. 17; 

Fig. 21 is a flowchart explanatory of a method by which the flow of 
process for evaluating the candidate vector and obtaining die direction of the noise 
magiietic field is used in the process for calculating the direction of the noise 
magnetic field in the flowchart shown in Fig. 17; 
10 Fig. 22 is a signal frequency spectrum diagram explanatory of how the 

operation of selecting a frequency for niaxijnizing a frequency spectrum is used in 
die process for calculating die direction of the noise magnetic field iii the flowchart 
shown in Fig. 13. 

Fig. 23 is a signal wavefomi diagram explanatory of the operation of 
15 specifying a period in which only a noise magnetic field exists through utilization 

of a fact diat the amplitude of the sensed magnetic field signal becomes small 

during die OFF period of die signal magnetic field that is turned OFF by a 

predetemiined procedure, in the calculation of the direction of the noise magnetic 

field according to die present invention mediod; 
20 Fig. 24 is a signal waveform diagram showing instantaneous variations in 

the amplitude of the sensed magnetic field signal when die signal magnetic field is 

periodically turned OFF in die calculation of the direction of the noise magnetic 

field according to die present invention method; 

Fig. 25 is a signal wavefomi diagram showing instantaneous variations in 
25 the amplitude of the sensed magnetic field signal when the signal magnetic field is 

randomly turned OFF in tiie calculation of the direction of the noise magnetic field 

according to the present invention mediod; 

Fig. 26 is a flowchart showing a method for calcidating the direction of the 

noise magnetic field according to Uie present invention in wliich the signal 
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magnetic field is turned OFF on tlie basis of a predetermined sequence, a period 
during whicli a particular magnetic field is OFF is specified by a sequence starting 
at tlie time when a correlation function between tiie sequence and tlie sensed 
magnetic field becomes maximum, and the direction of Hie sensed magnetic field 

5 in die specified period is regarded as tlie direction of the noise magnetic field; 

Fig. 27 is a flowchart showing a method for calculating the direction of tlie 
noise magnetic field according to tlie present invention in which the signal 
magnetic field is turned OFF on tlie basis of a predetermined sequence, and tlie 
starting time of a sequence indicating flie most likely ON/OFF state of the signal 

10 magnetic field is calculated from a plurality of times when a correlation fiinction 
between the predetermined sequence and tlie sensed magnetic field becomes 
maximum; 

Fig. 28 is a flowchart showing the metliod for calculating he direction of 
the noise magnetic field according to the present invention in which the signal 

15 magnetic field is turned OFF on tlie basis of a predetermined sequence, the 
correlation function between the predetemiined sequence and tlie sensed magnetic 
field is calculated at each of a plurality of time points at wliich Uie period of tlie 
sequence is equally divided, tlie magnetic field is projected at each time point to a 
vector tbmied by tlie calculated function, and one of vectors whose variance is 

20 minimum is regarded as the direction of the noise magnetic field; 

Fig. 29 is a flowchart showing another method for calculating the direction 
of the noise magnetic field according to the present invention in wliich tlie signal 
magnetic field is turned OFF on tlie basis of a predetermined sequence, a 
correlation function between the predetemiined sequence and tlie sensed magnetic 

25 field is calculated at each of a plurality of time points at which the period of the 
sequence is equally divided, tlie magnetic field is projected at each time point to a 
vector foniied by the calculated function, and one of vectors whose variance is 
minimum is regarded as the direction of the noise magnetic field; 

Fig. 30 is a flowchart showing tlie method for calculating tlie direction of 




the noise magnetic field according to the present invention in which the signal 
magnetic field is turned OFF on the basis of a predetermined sequence, and the 
starting time of a sequence indicating tlie most likely ON/OFF state of tlie signal 
magnetic field is calculated fi-om a plurality of time points when the correlation 
5 function between die predetermined sequence and the sensed magnetic field 
becomes maximum; 

Fig. 31 is a flowchart showing still another method for calculating the 
direction of the noise magnetic field according to the present iiwention in which 
the signal magnetic field is turned OFF on the basis of a predetermined sequence, a 

10 correlation function between tlie predetermined sequence and the sensed magnetic 
field is calculated at each of a plurality of time points into which the period of the 
sequence is equally divided, the magnetic field is projected at each time point to a 
vector fomied by the calculated function, and one of vectors whose variance is 
minimum is regarded as tlie direction of the noise magnetic field; and 

15 Fig. 32 is a perspective view depicting an example of a magnetic field 

sensing frame for use in tlie present invention. 

DETAILED DESCRIPTION OF THE INVENTION 
As depicted in Fig. 1, in a case where the position of a digging head 2 
under the ground surface 1, which is a signal magnetic field source to be sensed, is 

20 determined using a magnetic sensor 4 placed on the ground surface 1 at a proper 
position when a power line or similar magnetic noise source 3, which generates a 
noise magnetic field, is placed near the digging position .to be determined, there are 
present a signal magnetic field vector Hs provided fi^om the digging head 2 and a 
noise magnetic field vector Hn fi*om the magnetic noise source 3, such as power 

25 line. In this case, the magnetic sensor 4 senses a vector Hm that is a combined 
version of the signal magnetic field vector and the noise magnetic field vector 

Hn. 

Now, let the noise magnetic field of a position vector r and at time t be 
identified as a vector Hn(r, t). On the other hand, let the signal magnetic field 



generated by magnetic filed generating means for position sensing be identified as 
a vector Hs(r-rc, 0c, i)- Here, the vector 0c is an angle of orientation of the 
magnetic field generating means, which is defined by three angles of rotation in tlie 
coordinate system fixed to the ground that is the coordinate system fixed to the 
5 magnetic field generating means. 

Since die noise magnetic field and the signal magnetic field are botli 
sensed simultaneously and since the noise magnetic field vector Hn(t, t) varies 
randomly witli time, it is impossible to extract only the signal magnetic field vector 
Hs(r-rc, 0c, t) Irom tlie measured magnetic field vector Hm(r-re, 0c, t) unless 

10 statistical properties of tlie noise magnetic field are known and the noise magnetic 
field is signal-wise orthogonal to the signal magnetic field. Even if die statistical 
properties for separating the noise magnetic field fi^om the signal magnetic field are 
known prior to the deteraiiiiation of the position of die latter, the separation calls 
for a large amount of data, and hence tlie conventional scheme is not ever practical. 

15 According to the present invention, the direction of a vector en(r) of the noise 
magnetic field vector Hn(r, t) is obtained by separated means, and in a coordinate 
system shown in Fig. 2, a component vector HmV-rc, 0c, 0, projected to a plane 
vertical to die direction of a vector e„(r) of the measured magnetic field vector 
Hm(r-rc, 0c, t) as shown in Fig. 4 (S 1 , S2, S3). 

20 H./(r-r,, /)=H„Xr-r,, ^, /)- (H„,(r- r,, ^, O e„(r)X(r) (1) 
This component does not contain the noise magnetic field for die reason given 
below. Since 

H„/(r-.;, e^. /)=HXr-r,, O^^ /)+H„(r, /)= H^r- r,, 6>, 0+|H„(r, /)|e„(r). (2) 
it follows that 

25 H/(r-r., ^„ 0= H^r- r., ^„ i)-^i^-r^. />e„(r)>„(r> (3) 
fi^om which it can be seen that die projective component vector HmV^-rc, 0c, t) does 
not contain the component of the noise magnetic field vector Hn(r, t). 

However, die projective component vector Hm^(r-rc, 0c, t) has lost 
information of one axis by die projection on a plane vertical to the direction of the 

9 



vector en(r). That is, since the sanie projective components are obtained 
irrespective of tlie magnitude of a component parallel to the vector en(r), two 
independent components are obtained. 

Altliough any method can be used to obtain the two independent 
5 components; it is possible to use such a method as described below. 

Tliat one of coordinate axes of a measurement coordinate system Cm 
(which will be described later on), which is not parallel to the direction e„(r) of die 
noise magnetic field vector H„(r, t), is chosen. Let a unit vector in die direction of 
the chosen coordinate axis by identified as a vector e^. A vector product, 
10 ep^i=e,„xen(r), of die unit vector and the direction en(r) is perpendicular to die 
direction of the vector e„(r), and hence it is contained in the plane of projection and 
is perpendicular to the coordinate axis Cm. Let the magnitude of a vector obtained 
by projecting the projective component vector H,/(r-rc, Oc, t) in the direction of the 
vector Cpj, including the direction of the vector, be represented (S4) by a value of 
15 Hm,i V-rc, Gc, t). Tliatis, 

^nu^--;. ^c.O=H./(r-r,, ^i.O-e,.- (4) 
Cp, I = X e„(r> (5) 
The next step is to calculate a vector ep,2 perpendicular to die directions of 
the vectors Cp^ and en(r). Tlie direction of the vector Cp^ is also perpendicular to 
20 the direction of the vector ^(r), and hence it is contained in die plane of projection 
and is perpendicular to die direction of the vector Cp i as well. Letting the 
projection of the projective component vector Hm^(r-rc, 0c, t) in this direction be 
represented by Hn^^(r-rc, 0c, t), values H^,! V-fc, 0c, t)ep,i and Hm,2^(r-rc, 0c, t)ep;j 
are two independent vectors into which the projective component vector Hm^(r-rc, 
25 0c, t) is separated Here, 

^„.2'(r-r,, 0^, 0=HjXr-r,, 0., () e^:,, (6) 
ep,2 =ep^, xe„(r) (7) 

Then, by setting the position vector rc and angle-of-orientation 0c of the 
magnetic field source so tliat a theoretically calculated magnetic field, He(r— r©, 0c, t) 

10 



generated by tlie magnetic field source of the position vector r_substantially 
matches witli a projective component He^(r-rc, 0c, t) on tlie same plane as that of 
the measured magnetic field, it is possible to detect the position and orientation of 
the magnetic field source. 
5 The above description has been given of a case where the number of noise 

magnetic fields is virtually one, but when tlie number of noise magnetic fields is 
two, letting tlie direction of two noise magnetic fields Hnj(r, t), where j = 1, 2, be 
represented by vectors e„i(r) and en2(r), use is made of the projection of tlie 
measured magnetic field on a direction vector eN(r)=eni(r)x en2(r) in the coordinate 
10 system of Fig. 3 as depicted in Fig. 5. That is, 

H/(r^i;, 0^^ />H,,(r-r., ^, /)-(H„Xr-r,, ^, /) e^(r)>N(r) (8) 

is calculated (S 1 1 , S 1 2, S 1 3). hi the above 

e^(r)=e„.(r)xe„,(r) (9) 

where the symbol "x" indicates a vector product. Li this instance, an independent 
15 component iii tlie projective component is one, that is, only the magnitude of the 

vector. 

Now, let the coordinates of the magnetic field source be represented by a 
vector rc(x, y, z) and its angle of orientation by a vector 0c(O>t, 0y, 9z)- In the 
following description, 9x, Gy and will be referred to as an angle of rotation, a tilt 
20 angle and an azimutli, respectively. When the signal magnetic field is axially 
symmetric, tlie axis of symmetry is regarded as the x- axis, and the angle of 
orientation is set to a vector 9c(9y, G^^. 

Further, according to tilie present invention, it is possible to calculate the 
magnitude Hni(r-rc, 9c) of each component of the original signal magnetic field 
25 vector Hni(r-rc, 9c, t) by synchronous detection of the measured magnetic field 
vector Hni(r-rc, 9c, t) tlirough die use of a proper one component in Eq. (1) or (6) 
that is a projective component. 

HXr-n, ^e) = (H„Xr-r,, ^,,/K>-r,, z))^. (10) 

In the above, Hn,V"rc, 9c, t) is a proper component of the projective 
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component given in Eq. (1) or (6). 

Fig. 4 shows die case where tlie number of noise magnetic field sources is 
essentially one. The processing of Fig. 4 is repeated at each measuring point. In 
the case of essentially two noise magnetic fields, too, the processing shown in Fig. 
5 5 is repeated at each measuring point as is the case with Fig. 4. 

The magnitude Hn,(r-rc, ©c) of tlie original signal magnetic field Hni(r-rc, 0c, 
t) can be obtained by calculating Eq. (10) after calculating the projective 
component at each measuring point. Concretely, tlie processing depicted in Fig. 
6 or 7 is carried out [(SI 1, S12, S14, S15, S16) or (Sll, S12, SMa, S16)]. 
10 How to calculate the position and angle of orientation: 

(1) Definitions of tlie coordinate system and the angle of orientation: 

A definition will be given first of tlie coordinate system necessary for the 
description of the invention. 

A coordinate system is set which is fixed to the earth with die z- axis in a 
15 vertical direction (upward), wliich system will hereinafter be referred to as a 
measuring coordinate system vector Cm- Tlie x- and y- axes are properly set so 
tliat they forni a right-hand system. For example, they are set in parallel to die 
direction in which die side of a measuring fi*anie is projected on a horizontal plane. 

This has for its object to calculate the coordinate vector re(x, y, z) and 
20 aiigle-of orientation vector 0c(6x, 9y, 9z) of the magnetic field source in die 
coordinate system. 

On the other hand, as the coordinate system vector Q of the magnetic field 
source, a coordinate system with die axial direction of the magnetic field source as 
die X axis is set so that die y- and z- axes are horizontal and vertical (upward), 
25 respectively, when the magnetic field source is placed horizontally. 

The aiigle-of-orientation vector 9c of the magnetic field source is defibned as 
the angle of rotation between the measuring coordinate system vector Cm and die 
coordinate system vector Q as described below. In the first place, a coordinate 
system Cco parallel to the measuring coordinate system vector Cm is turned by an 
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aziniutli aiigle G^. about tlie z- axis (in eitlier one of the coordinate system). This 
coordinate system will be referred to as a coordinate system vector Cd. Next, the 
coordinate system vector Qi is tumed by a tilt angle 9y about the y- axis of the 
coordinate system vector Qi itself Tliis coordinate system will be referred to as 
a coordinate system vector Cc2. Further, the coordinate system vector Cd is 
turned by an angle of rotation 9z about the x- axis of the coordinate system vector 
Cc2 itself The angle-of-orientation vector 0c is determined so that the resulting 
coordinate system becomes the coordinate system vector Q. 

(2) Explanation of independent measurands and unknowns: 
When the number of noise magnetic fields is essentially one, the 
measuiement of tlie magnetic field at one place will provide two independent 
measurands. When tlie number of noise magnetic fields is essentially two, the 
measurement of the magnetic field at one place will provide one independent 
measurand. Further, in case of obtaining components of the original signal 
magnetic field by synchronous detection, the measurement of the magnetic field at 
one place will provide three independent measurands. On tlie other hand, three 
coordinate components of the coordinate vector rc(x, y, z) of the magnetic field 
source are unknowns. The azimuth angle Oj is not obtainable without a different 
criterion such as the direction of the earth magnetism, and hence it is also an 
unknown unless the digging head is provided with a direction sensor When the 
digging head is made of a magnetic material, or when a buried steel pipe or similar 
magnetic body lies near the digging head, an accurate direction cannot be obtained 
even if the direction sensor is provided; hence, the azimuth angle 6z is an unknown 
in many cases. The tilt angle 9y can easily be detected by a tilt angle sensor that 
detects the vertical direction, and hence it is known in many cases. The same is 
true of the angle of rotation Q^. In particular, when the signal magnetic field is 
axially symmetric, if tlie axis of symmetry is regarded as the x- axis, the angle of 
rotation becomes meaningless and hence can be ignored. 

At any rate, measured magnetic fields need only to be obtained at differctnt 
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positions so that independent measurands larger in number tlian unknowns. 
Arrangement of the measuring system: 

For example, as shown in Fig. 1, magnetic fields are measured using a 
required number of tliree-axis magnetic sensors disposed on the ground so tliat 
5 their relative positions are known. Since the direction of the noise magnetic field 
changes for each position of measurement, it is necessary to perfomi for each 
position the step (S21) for detecting the position of tlie noise magnetic field and the 
step (S22) for calculating die projective component of the measured magnetic field 
as shown in Fig. 8 that is a flowchart of tlie procedure for position determination 
10 by the present invention. In case of calculating tlie component of the original 
signal magnetic field by synchronous detection, the processing therefor (in a 
specified step) needs to be carried out for each position. 

Flow of measurement processing: 

A detailed description will be given of the detection (S21) of the direction 
15 of tihe noise magnetic field, the calculation (S23) of the position of the signal 
magnetic field source tlirough the use of the projective component and the 
calculation of tlie position f of the signal magnetic field througli the use of the 
signal magnetic field component. 

[EMBODJMENTS] 
20 Embodiment in the case of virtuallv one noise magnetic field: 

In case of using the projective component, when the number of unknowns 
is Nu(>l) and tlie number of noise magnetic fields is virtually one, the projective 
component vector Hn,V~^^c, ^c, t) expressed by Eq. (1) is calculated firom magnetic 
vectors Hni(rr~rc, Oc, t) measured at Nu/2 or more different positions, and the 
25 position vector rc and the angle-of-orientation vector 0c are detemiined (S31) as 
depicted ki Fig. 9 so tliat a projective component vector at each magnetic field 
sensing point is essentially equal to the projective component vector He'^Cr-rc, 9c, t) 
of a tlieoretically calculated magnetic field of the projective component at each 
position of measurement, by which it is possible to detect (S32) the position and 



orientation of the magnetic field source. 

In case of using the magnitude of tlie signal magnetic field synchronously 
detected by the projective component, tlie magnetic field is measured at Nu/3 or 
more diflferent positions. When the number Nu(^l) of unknowns is even, the 
5 number of unknowns and the number of independent measurands can be made 
equal to each other; if setting magnetic field vector Hni(r-rc, Gc, t) is measured at 
positions Nm=Nu/2, a Ny number of such equations as given below need only to be 
solved. 

K./(r.-r,,^.,OX-^e./(i-.-r,. ^c)=0>^ (H) 
10 where q = 1, 2 and represents two directions parallel to the plane of projection but 
parallel to each otlier (S33). Accordingly, Hni,q''(r-rc, 0c, t) and He,q^(r-rc, 9c, t), 
where q == 1, 2, are tlie magnitudes of q-direction components of the measured 
magnetic field vector Hmlr-rc, 9c, t) and the tlieoretical calculated magnetic field 
vector Hc(r-rc, 9c, t), respectively. Let it be assumed that the vector 9c represents 
1 5 any one of 9x, 9y, 9^, (9y, 92,), (9^, 9x), (9^, 9y), (9^, 9y, 92) and ^ and that the vector Vc 
represents any one of x, y, z, (y, z), (z, x), (x, y), (x, y, z) and (|), where ^ represents 
an empty set 

For example, when unknowns are the position vector rc(x, y, z) and 
azimuth angle 9z of the magnetic field source, the magnetic field is to be measured 

20 at two different positions and four equations such as given below are solved, by 
which it is possible to obtain (S34) the position vector rc(x, y, z) and azimuth angle 
9z of the magnetic tield source. 

K;0\-re, ^z>O),-^./0v-';. 0 = 0, = 2;^=1, 2. (12) 
hi the above, < . >t represents a time average. 

25 As depicted in Fig. 10, when the magnetic field is measured at Nm (>Nu/2) 

different places more than Nu/2 in excess of tiie number Ny of unknowns, since a 
larger number of independent measurands than the number of unknowns can be 
obtained (S41, S42, S43), the position rc(e, y, z) and the angle of orientation 9z are 
calculated wliich provide 
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where < . >t represents tlie time average and min { .} means that vectors rc and 0c 

are changed to obtain vectors rc and 9c tliat minimize tlie contents in { . }. Further, 
a symbol Wk,q indicates weighting. Hie Eq. (13) can be replaced witli the 
5 following equations (S44). 

,|^(k™./(»lc-re. 0.,jt)^ - //, /(r, - r,, (14) 
min|£i:>v,,J(|//^;Xr,-.;, O^, ^H) ^ - ^ (r, - v^, ffji (15) 



Assume that vectors rc and 0c in Eqs. (13), (14), (15) and (16) have the same 
10 meaning as in the case of Eq. (11). 

hi the above description it does not matter whether the measured magnetic 
field vector Hn,(r-rc, 0c, t) is a signal having passed through a band pass fQter that 
permits tlie passage tlieretlirougli of only components close to the frequency of,the 
signal magnetic field or a wide-band signal that is inhibited from the passage 
15 through the band pass filter, but the use of the signal having passed through the 
band pass filter increases the possibility of determining the position of the 
magnetic field with high reliability. 

(Embodiment in the case of essentially two noise magnetic fields) 
When the number, of unknovms is Nu (^1) and the number of noise 
20 magnetic fields is virtually two, tlie projective component vector HmV-rc ©c, 0 
expressed by Eq. (4) is calculated from magnetic vectors Hn,(r-rc, 0c, t) measured 
at Nu/2 or more different positions, and the position vector rc and the 
angie-of-orientation vector 0c are determined (S51, S52) as depicted in Fig. 11 so 
that the above-mentioned projective component vector obtained at each position of 
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measurement esseatially matches witli the projective component vector He''(r-rc, 9c, 
t) of a tlieoretically calculated magnetic field at each position of measurement, by 
which it is possible to detect (S53) the position and orientation of the magnetic 
field source. 

hi this case, the number Nu (>1) of unknowns and the number of 
independent measurands can be made to be equal to each otlier witliout fail; the 
magnetic field vector Hm(r-rc, Gc, t) is measured at different positions of the same 
numbpr as that Nu number of unknowns, and Nu number of such equations given 
below need only to be solved. 

{Hj(r^, - r„ 0,, t))^ - Hfiry - i; , 0j= 0, k = l ... (17) 

In this case, -symbols H,„V-rc, %, t) and hJ(t-To, 0c, t) are the respective 
magnitudes of die projective component vector H„/(r-rc, 9c, t) and the vector 
HeV-fc, ©c, t) of the measured magnetic field vector Hm(r-rc, 9c t) and the 
theoretical calculated magnetic field vector He(r-rc, Oc, t). Let it be assumed that 
the vector 0e represents any one of 9,, 9y, 9^ (9y, 9,), (9„ 9,), (9,, 9y), (9^, 9y, 9^) 
and ^ and tliat the vector represents any one of x, y, z, (y, z), (z, x), (x, y), (x, y, z) 
and (|), where ^ represents an empty set. 

For example, when unknovsms are the position vector rc(x, y, z), azimuth 
angle 9^ and tilt angle 9y of tlie magnetic field source, the magnetic field is to be 
measured at five different positions and four such equations given below are 
solved, by which it is possible to obtain tlie position vector rc(x, y, z), azimutli 
angle and tilt angle 9y of the magnetic field source. 

{f^J'(ry-r.,0o,d),-ff/(r^-r,,0j=O,k = l,...,5. (18) 
In the above, < . >t represents a time average. And the vector 9c = 9c(9y, Q^. 

As depicted in Fig. 12, when the magnetic field is measured at Nm (>Nu) 
different places more tlian Nu in excess of the number Nu of unknowns, since a 
larger number of independent measurands than tlie number of uiiknovms can be 
obtained (S61, S62), the position rc(e, y, z) and the angle of orientation 92 are 
calculated which provide 
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(19) 



where < . >i represents the time average and min { .} means that vectors and 0^ 



are changed to obtain rc and 0c tliat minimize the contents in { . }. Further, a 
symbol Wk ,, indicates weighting. The Eq. (1 9) can be replaced with tlie following 
5 equations (S63). 



Assume that rc and 0c in Eqs. (19), (20), (21) and (22) have the same meaning as in 

10 the case of Eq. (16). 

In die above description, it does not matter whether the measured magnetic 
field vector H„,(r-rc, 0c, t) is a signal having passed through a band pass filter that 
permits the passage theredirough of only components close to the fi-equency of the 
signal magnetic field, or a wide-band signal that is inhibited fi^om the passage 

15 through the band pass filter, but the use of the signal having passed through the 
band pass filter increases the possibility of determining the position of the 
magnetic field with high reliability. 

Further, in case of using the signal magnetic field component obtained by 
synchronous detection, the direction of the noise magnetic field is determined and 

20 is used to obtain the projective magnetic field; and the subsequent processing is 
common to die cases of one and two noise magnetic fields. Hie projective 
component vectors H„,''(r-rc, 0c, t) at respective places are calculated fi^om 
magnetic field vectors Hn,(r-rc, 0c, t) measured at Nu/3 or more different places, 
and a proper one of tlie calculated vectors is used as a reference signal to perform 




(20) 




(21) 




(22) 
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synchronous detection of the measured magnetic field vectors Hm(r-rc, 9c, 0, 
thereby obtaining the magnitude Hs(r-rc, Gc) of the original signal magnetic field 
component. By detennining tlie position vector rc and tlie angle-of-orientation 
vector 0c of the ma^etic field source so that the magnitude of the original signal 
5 magnetic field component and the magnitude Hc(r-rc, Gc) of the theoretical signal 
component are equal to each other, it is possible to detect tlie position and 
orientation of tlie magnetic field source. 

When tlie number Nu (>I) of unknowois is a multiple of 3, tlie number of 
unknowns £ind tlie number of independent measurands can be made to be equal to 

10 each other; if tlie magnetic field vector Hm(r-rc, Gc, t) is measured at positions 
Nm=Nn/NV3, Nu number of such equations as given below need only to be solved. 

HXr. - ^c)- H^r. - 0= 0, * = 1, ... (23) 
Fig. 13 shows tlie flow of processing (S61, S62, S65). Here, Hs(r-rc, Gc) and 
Hc(r-rc, Gc, t) are tlie signal magnetic field calculated fipm the measured magnetic 

15 field vectors Hm(r-rc, Gc, t) and tlie magnitude of the theoretical calculated 
magnetic field vector, respectively. The magnitude Hs(r-rc, Gc) of the signal 
magnetic field vector is an averaged quantity as already explained. Let it be 
assumed that the vector Gc represents any one of Gx, Gy, G^, (Gy, G^, (G^, G^), (G^, Gy), 
(Gx, Gy, Gz) and ^ and that the vector rc represents any one of x, y, z, (y, z), (z, x), (x, 

20 y), (x, y, z) and <}), where <|) represents an empty set. 

For example, when an unknown is the position vector rc(x, y, z) of the 
magnetic field source, the magnetic field is to be measured at one place and three 
such equations as given below are solved, by which it is possible to obtain the 
position vector rc(x, y, z) and azimuth angle Gz of the magnetic field source. 

25 H,(r-rJ-HXr-rJ=0. (24) 

Tlie vector r is the vector of tlie position of measurement. 

When the magnetic field is measured at Nr, (>Nu/2) different places more 
than Nu in excess of the number Nu of unknowns, since a larger number of 
independent measurands than tlie number of unknowns can be obtained, tlie 
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position rc(e, y, z) aiid tlie angle of orientation 0^ are calculated which provide 
'^ij]|Z^>,Ol - r., 0^- H,(r, - r , , (25) 

where min { .} means that vectors and 0c ai e changed to obtain vectors rc and 0c 

that minimize tlie contents in { . }. Furtlier, a symbol Wk,q indicates weighting. 
5 The Eq. (25) can be replaced with the following equations. 

min{Ev-*.(lHsOlc - re, ^cj- K(r. - ^.]||)'|. (26) 

™ii|{£^jHXr. - r., 0^- H^r, - r„. ^J'j. (27) 

Fig. 14 shows flie flow of processing (S61, S64, S66). 

A description will be given below of how to determine the direction of the 
10 noise magnetic field in the two embodiments described above. 

(How to detemiine the direction of the noise magnetic field in the case of 
virtually one noise magnetic field) 
(First metliod) 

The first method for detemiining the direction of the noise magnetic field is 
15 a method where, in case of no signal magnetic field, a noise magnetic field is 
measured through the use of the same measuring system as that for measuring the 
signal magnetic field. This situation is such as shown in Fig. 15 in which the 
detection of tlie digging position is disturbed by a noise magnetic field source lying 
in the vicinity of tlie digging route. Another case is that where the signal 
20 magnetic field source is equipped with a function of receiving a commmd sent, for 
example, from the ground by some means and responsive to the command to stop 
the generation of the signal magnetic field. 

In tliis instance, as depicted in Fig. 16, letting the measured magnetic field 
be represented by a vector Hm(r-rc, 9c, t), since this is essentially a noise magnetic 
25 field vector Hn(r, t) as shown at a step an average value of its absolute values 
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can be used (S72) to calculate a direction vector, e„(r)=en,x(r), e„,y(r), e„,7.(r)), of tlie 
noise magnetic field as follows: 

^",a('-)=-r==T= = ^*==y, «=^,j.,z. (28) 

Alternatively, a root-mean-square value of the above absolute values can be used to 
5 determine the direction of tlie noise magnetic field by 

In this case, a symbol Hm,a(r-rc, Oc, t) is an a component of the measured magnetic 
field (die noise magnetic field), and a is any one of x, y and z. 

In the above description it does not matter whether the measured magnetic 
10 field vector Hm(r-rc, 0c, t) is a signal having passed through a band pass filter that 
permits the passage therethrough of only compon^ts close to the fi-equency of tlie 
signal magnetic field or a wide-band signal that is inhibited from the passage 
through tlie band pass filter, but the use of the signal having passed through the 
band pass filter increases the possibility of determining the position of the 
15 magnetic field with high reliability. 

(Second method) 

Step 1 

As shown in Fig. 17, at first , the frequency spectrum Hn,((») of the 
measured magnetic field vector Hm(r-rc, Gc, t) is calculated (S81, S82, S83) by the 
20 following equation. 

H = F (H„ (r, - r., a, 0) (30) 

In this case, a symbol F( . ) represents a Fourier transform; the three components x, 
y and z of die measured magnetic field vector H„,(r-rc, Gc, t) are each Fourier 
transformed. In practice, the above-mentioned fi-equency spectrum can be 
25 calculated by FFT (fast Fourier transform) or the like of sampled values of the 
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measured magnetic field vector Hni(r--rc, Gc, t). 
Step 2 

Next, an angular frequency coj, where i = 1, 2, Ns, of a large-amplitude 
component, such as a line spectrum, is selected (S84) firom the absolute value 
5 |Hin((o)| of die jfrequency spectrum. For the component of each angular frequency 
(Di, where i = 1, 2, . . . , n,is, a candidate unit vector en(r, (Dj), where i = 1, 2, . . , Nns, of 
the direction of die noise magnetic field, is calculated (S85) by die mefliod (1) or 
(2) described below. 

(1) The absolute values of the Fourier-transformed x, y, and z components 
10 of die angular frequency concerned are used to calculate the candidate unit vector 

en(r, coi), where i = 1, 2, . . . , Nns, of die direction of the noise magnetic field, by the 
following procedures: 

e„lr,6>J= 1^1^ . ^2 ' / = 1, ... , TV... (31) 

\ a-x. y, z 

where H,„a((Di), where a = x, y, z and i = 1, 2, N„s, is a (Di component by the 
1 5 Fourier transform of the a component of the measured magnetic field. 

(2) A narrow-band filter is formed whose pass band uses, as the center 
fi-equency, the angular fi^equency ©i, where i = 1, 2, n„s, and the same method 
as by Eq. (21) or (22) is used to calculate the candidate unit vector en(r, coO, where i 
= 1 , 2, . . . , N,K . That is, the candidate unit vector e„(r, (jOi)= (e„pc(r), en^(r), e„,2(r)), 

20 where i = 1 , 2, . . . , Nns ,is calculated by 



a = x,y,z;i = \ ... , N^. (32) 



or by 

, ^ J((^n..,('-k-';. ^c,^„ 01) 



(33) 
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step 3 

The candidate unit vector en(r, a)i)= (e,^(r), e,^y(r), e,^(r)), where i = 1, 2, . . . , 
N„s , is considered as tine direction vector e„(r) of the noise magnetic field, and for 
each angular frequency ©i, where i = 1, 2, . . . , Nns , the same method as by Eq. (1) 

5 is used to calculate a projective component vector Hra'*(r— rc, 9c, oOj, t). 
H/(r-r,, e^, H.„(r-r„, 0^, t) 

-(H...(r-r^, e^,t).^M ^,)>„(.% o>^^ = \, ... , Ks- (34) 
The reason for which the angular firequency ©i is contained as a variable of the 
projective component is to explicidy point out that the projective component is 

10 dependent on the angular fijequency coi, where i = 1, 2, • A roper time 

interval Ttest, which consists of Nust durations Tg^ where k = 1, 2, Hest, each 
having a short time length At, is chosen, and the variation of the projective 
component vector Hn,^(r-rc, 0c, co;, t) for each duration Ti^k, where k = 1, 2, . . . , Hest, 
is evaluated (S86). Assume, here, that each duration Tj^^, where k = 1, 2, . . . , Nusi 

1 5 does not overlap other durations. Concretely, a variance of Niest statistics Vevaijc(o>i) 
of the Ntest, where k = 1, . . . ^ Hesu which are calculated by any one of the methods 
described below, is calculated. 

(1) One or botli of tlie means of absolute values of two orthogonal components by 

v.™...W=K. 0,.,, t\^,g=h 2;k=l, ...,N,^; /=1, ... , AT,. (35) 

20 where {)^^ ^ represents the mean value in the duration Tg^ and Vevaijc is a statistic 
calculated for tlie duration Ti^^. 

(2) Mean of absolute values 

v^^,, i,(co} = ^U J (r-r^, ^, a>„ i^^ ^ k = l ... , N^; /= 1, ... , N„^. (36) 

(3) One or both of the means of squares of two orthogonal components by 
25 ,(a,,)=((/f„./(r-r,. 0^, o>,, i)^) , q = l, 2- k = h ... , N^; i=l, ... , N^. 

(37) 
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(4) One or both of square roots of tlie means of squares of two orthogonal 
components by 



(^i) = -^((^n„/('-rc. e., , q = l,2-k = I, ... , N,^; / = 1. ... , 



(38) 

5 For tlie statistics Vevai,k, where k = 1 , . . . , Niest, calculated by tliese equations, the 
following equation 

Jmeank(j[v^3i, y,(co, )- mean,(v,,„, ^(o?^ Vif^ 
varW=JI ^ J (39) 

is calculated to obtain (S86) a value of (Di,,„in that is die angular frequency which 
minimizes var(coi). hi the above, meank(.) indicates averaging for the suffix k, 
10 that is, 

mean,(.) = ^-^. (40) 

The magnetic field of tlie angular frequency (Dynin derives from tlie noise 
magnetic field, and tlie direction of tlie noise magnetic field becomes a vector ^,(r, 

15 Incidentally, the angular frequency coi^miib which minimizes var((ji)i), needs 

only to be measured at one place and need not be obtained at every place where to 
measure tlie magnetic field. 

Witli tliis metliod, it is also possible to calculate a fluctuation in the 
direction of a vector Hm**(r-rc, 9c, coi, t) as well as an amphtude fluctuation given by 
20 Eq. (39) and to select die angular frequency cOj^min at which the direction fluctuation 
becomes minimum or smaller dian a predetermined value. 

hicidentally, in the above description, die measured magnetic field vector 
Hni(r-rc, Gc, t) in Step 1 is a wide-band signal, and in Step 3 it does not matter 
whetlier the measured magnetic field vector Hm(r"rc, ©c, t) is a signal having passed 
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tlirough a band pass filter that permits the passage therethrough of only 
components close to tlie frequency of the signal magnetic field or a wide-band 
signal that is inliibited from tlie passage tlirough tlie band pass filter, but tlie use of 
the signal having passed through the band pass filter increases the possibility of 
5 detemiining the position of the magnetic field with high reliability. 

Fig. 18 shows how to select frequencies fi(=CD 1/271), f2(=o)2/27i), 
fii(=(Q)n/27i, where n=Nns), at which the frequency spectrum becomes maximum; 
Fig. 19 shows the flow of processing including step S91 for obtaining the 
candidate vector; and Figs. 20 and 21 show tlie flow of processing including steps 

10 SlOl and S102 or steps Sill and S112 for evaluating the candidate vector and for 
detecting the direction of the noise magnetic field. 

The fiequency spectrum Hni(o) need not always be used. That is, tlie 
signal magnetic field is periodically turned OFF/ON following a predetenniiied 
procedure; the period Tpenod is divided into equally-spaced durations; the candidate 

15 unit vector e„(r, tj), where i = 1, . . . , Nns, is used in place of the candidate unit vector 
en(r, oOi), where i = 1, Nns, which is calculated by Eqs. (25), (26) and (27); and 
tliereafter, the duration tliat minimized the variance by Eq. (33) is calculated by the 
processing described above. By tliis, the vector en(r, tj) in that duration can be 
adopted as tlie direction of he noise magnetic field. 

20 (Third metliod) 

In the second method, the large-amplitude angular frequency coi, where i = 
1,2, . . . , Ns, is selected from the absolute value |Hni.(cD)| of the frequency spectrum 
Hm(Q)), but the vector e„(r, cos^in) can be obtained as the direction of the noise 
magnetic field in exactly the same mamier as in tlie case of the second method, by 

25 selecting a proper frequency band neighboring the frequency of die signal 
magnetic field, setting properly-spaced test frequencies free from the frequency of 
the signal magnetic field in the frequency band and regarding the test frequencies 
as the angular frequency (0\ in the second method. 

As is tlie case with tlie second method, the angular frequency at 
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which var(a)i) becomes minimum needs only to be obtained. 

Fig. 22 shows how to select die frequencies fi(=a)i/27cX f2(=C02/27r), 
fn(=(a)n/27i:, where n=N,^) at which die frequency spectrum becomes maximum. 
The flow of the subsequent processing is the same as depicted in Figs. 17, 19 and 
5 20. 

(Fourtli method) 

Tlie signal magnetic field is periodically stopped under control of a 
predetennined procedure. For example, tlie signal magnetic field is periodically 
stopped by a predetemiined time interval. Since the intensity of the magnetic 

10 field being measured decreases wliile tlie signal magnetic field is stopped, the OFF 
period of the signal magnetic field is identified by regarding the 
intensity-decreasing period essentially as tlie predetemiined OFF period, and tlie 
direction of the magnetic field measured during the OFF period is used as the 
direction of the noise magnetic field. The direction of the noise magnetic field 

15 can be obtained using tlie same method as the first one. Fig. 23 shows how. die 
amplitude of die measured magnetic field in this method varies with time. 
(Fiflfa method) 

The signal magnetic field is periodically stopped following a 
predetermined procedure. This is carried out as described below under (I) and 
20 (2). 

(1) To stop the signal magnetic field on a rectangular-wave-wise: / 

Tlie signal magnetic field is repeatedly turned ON and OFF widi the period 
Tperiod, for instance. 

25 and when a sequence s(t) is 1, the signal magnetic field is tjmed OFF, but when 
the sequence is -1 , die signal magnetic field is turned ON. In the above, 

^kM ~ = Vrk>d5 ^ = U 2, 3, ... (42) 

(2) To stop die signal magnetic field on a pseudo-random-signal-wise 
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basis. 

For example, when tlie value is "-1" in a random sequence like an 
M-sequence consisting of unit periods T„„it of llie same lengUi Nm, the signal 
magnetic field is tumed ON, but when the value is "1," the signal magnetic field is 
5 tumed OFF; and this sequence is repeated. In this case, the time average of the 

sequence is set to 0. 

Fig. 24 shows temporal variations in the amplitude of flie measured 
magnetic field in the case of the signal magnetic field being stopped by the metliod 
(1). Fig. 25 shows temporal variations in the amplitude of the measured magnetic 
10 field in the case of the signal magnetic field being stopped by the method (2). 

Next, as depicted in Fig. 26, tlie correlation function between the sequence 
s(t) and the nomi of tlie measured magnetic field or the absolute value of its 
particular component is calculated (S121, S122, S123). As the correlation 
function, it is possible to use any one of those calculated by 

15 Rir)=C;''''nKir-r^, ^^^^ 

j;^^'^^^-^(H„(r-r.. ^, t)ysii-r)d, (44) 
/^(r)=f;f''^''*^K„(r-re, 6»,/K'-rV/,a = x,3;,z. (45) 

7^(r)=f ^"^'-"J(//„,„(r-r,, 0^, Ol a = x,y,z. (46) 

hi this case, tlie period for which to detect the correlation is set to an integral 

20 multiple NrTperiod of die period Tperiod- 

Tlie OFF state of the signal magnetic field can be detected (S125) from the 
time T=tsy„c at which any one of the above conelation fijiictions become maximum 
(S124). That is, a sequence s(t-tsy«:), wliich starts at the time x=tsync, is used, and 
when the sequence s(t) is "I," die signal magnetic field is regarded as being OFF; 
25 in this way, the ON/OFF operation of tlie signal magnetic field is determined. 
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Ill tiie tlius detemiined signal magnetic field OFF period the direction of 
the measured magnetic field is dfeteqted, and flie direction is regarded as the 
direction of tlie noise magnetic field (S126). Wlien tlie sequence s(t) is (a 
second numerical value), the signal magnetic field is turned ON, and when the 
5 sequence is "1" (a first numerical value), the signal magnetic field is turned OFF; 
this sequence is repeated in this way. 
(Sixth method) 

As is the case witli the fifth method, the signal magnetic field is 
periodically stopped under control of such a predetermined procedure as described 
10 below. 

(1) To stop the signal magnetic field on a rectangular-wave-wise: 
The signal magnetic field is repeatedly turned ON and OFF with tiie 
period Tperi«xb for instance. 

5(0 = 1, /*</</t+W (41) 



= -1, tk + t^op^t^(M- 



15 and when a sequence s(t) is 1, the signal magnetic field is turned OFF, but when 
the sequence is -1, the signal magnetic field is turned ON. In the above, 

(^..-t, = T^^,k=\,2,3,... (42) 

(2) To stop tlie signal magnetic field on a pseudo-random-signal-wise 

basis: 

20 For example, when the value is "-l" in a random sequence like an 

M-sequence consisting of unit periods Tunit of the same length Nm, the signal 
magnetic field is tumed ON, but when the value is "1," the signal magnetic field is 
tumed OFF; and this sequence is repeated. 

In this case, the sequence s(t) is chosen so that it changes for each 
25 predetermined time unit Atu„ii. And, the time average of tlie sequence is "0." 

Next, as depicted in Fig. 27, tlie correlation function between the sequence 
s(t) and the norm of the measured magnetic field or the absolute value of its 
particular component is calculated (S131. S132, S133) as in case of flie fifth 
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method. As tlie correlation fiuiction, it is possible to use any one of those 



calculated by 

;e(,). J'J-^^-|h (r_ /)|K'-^>/'. (43) 

Rir) = J;""^'--VK (■•-••«' ^' 0)> - (44) 

5 , ^Ur)=j;;''''''^K„(r-n, 6>,4</-rK« = w. (45) 



In this case, the period for which to detect the coirelation is set to an integral 
multiple NxTperiod of tlie period Tpcdod- 

In this instance, tliere are present, in general, plural times T=%y^^ (k = 1, 2, 
10 . . . , Nsync) in which the correlation function becomes maximum and the maximum 
value exceeds a predetermined value (S 134). Assume, for example, that tsyncjt, (k 
= 1, 2, Nsync) is an arrangement of such times in order of time. When the 
correlation value between the sequence s(t) and the signal magnetic field is 
appropriate, 

15 A/^„^ k = ',y«^ k - '»yno. 1 , k=2, .... N^„^. (47) 

is virtually an integral multiple of the time unit At^it. Thai, the average of the 
value resulting fi-om the subtraction of an integral multiple Msync^cAtunit of the time 
unit AUync, where k = 2, Nsy™:.k, from At^;^, where k = 2, Nsy„c, is 
calculated by 

X (A^syno, k -M=y„c, k^unu) 

20 ^.y..=-^ Tj— . (48) 

In this case, 

provides the beginning of the sequence signal corresponding to the ON/OFF 
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operation of the signal magnetic field. 

Accordingly, the period during which the signal magnetic field is OFF can 
easily be set based on the sequence s(t-tsync). 

By applying the same method as die first method to the magnetic field 
5 vector Hm(r-rc, 0c, t) measured in this period, the direction e„(r) of the noise 
magnetic field vector Hn(r, t) can be calculated (S 1 38). 

(Seventli method) 

This mediod will be described below with reference to Figs. 28, 29, 30 and 

31. 

10 As is the case with the fifth method, the signal magnetic field is 

periodically stopped, for example, by such a procedure as described below. 

(1) To stop tlie signal magnetic field on a rectangular-wave-wise: 

The signal magnetic field is repeatedly turned ON and OFF with the 
period Tpenod, for instance. 

a (0 = 1, 4 < / < 4 + / . . 
15 V/ ' * A „„p 

= -1. '*+'s.op^'^'*.,- . ^ ^ 

and when a sequence s(t) is 1, the signal magnetic field is turned OFF, but when 
the sequence is -1, the signal magnetic field is tumed ON. In the above, 
'k+i-4 = ?;=riod, ^=1, 2, 3, ... (42) 

(2) To stop Hie signal magnetic field on a pseudo-random-signal-wise 

20 basis: 

For example, when the value is "-1" in a random sequence like an 
M-sequence consisting of unit periods Tu^i of tlie same length Nm, the signal 
magnetic field is tumed ON, but when the value is "1," the signal magnetic field is 
tumed OFF; and this sequence is repeated accordingly. In this case, the time 
25 average of tlie sequence is "0." 

Next, the correlation fimction between the sequence s(t) and the measured 
magnetic field Hm(r-rc, Oc, t) is calculated (S141, S142, S143). The period Tpenod 
is divided into equally spaced Ndiv sections of a length Tdiv, and either one of the 
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following calculations is conducted. 

C air - - k = 1, ... . ;V,,; a = x,y,z. (50) 

In tliis case, a symbol R<x(tk), where a = x, y, z and k = 1, Ndiv, is the time 
5 correlation between an a component Hm,a(r-rc, 6c, t) of tlie measured magnetic 
field Hm(r-rc, 0c, t) and tlie sequence s(t). And 

^-7',,,, k = \...,N^,. (52) 
The measured magnetic field may also be correlated with die time that is an m- 
multiple of the period Tpenod- That is, 

10 /^(rJ=J,^"""'"k.„(r-r„, 6»,/K0-'.^V> ^=^. ••■.^*v; ct = x,y,z. (53) 

KiO=\':J""'^iQi^Xr-r., 0.. OyS^^(t-t,)dt, k = l, ...,N,„ a=x,y,z. (54) 

In tliis instance, the sequence s(t) is replaced with Smp(t), where 
5.p(0-s(/ + 7;^^) ^^^^ 
5.„p(/) = .s<0, o<t<7;,,,^. 

This is followed by calculating the component Hm**(r-rc, 9c, tk, t), where k = 
15 1 , . . . , Ndiv , of the measured magnetic field Hm(r-rc, 9c, t) projected on the vector 
e„(tk)=(Rx(tk), Ry(tk), Rz(tk)), where k = 1, . . . , Ndiv 
formed by correlation functions Ra(tk), where a = x, y, z , corresponding to the 
respective components x, y and z of the measured magnetic field Hn,(r— rc, 9c, t). 
Here, a symbol tic contained as a variable of tlie projective component indicates that 
29 the projective component depends on the variable t^. 

hi tliis metliod, the vector e,i(tk) can be used as the direction of the noise 
magnetic field through utilization of tlie time tk in which a fluctuation in the 
absolute value of the projective component vector Hm(r-rc, 9c, tk, t), where k = 1, 
. . . , Ndiv, 
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j((H„,'(r-re, 0^, /u. /)-(|HJXr-re, B., t^, /]|>J) 

becomes minimum or smaller tlian a predetermined value (SI 45a). In liie above 
means tlie calculation of the time average. 

5. Further, variances of the x-component Hn,,x(r-rc, Oc, tk, t), y-component 

Hnij,(r-rc, 0c, ^, t) and z-component H„^(r-rc 8c, tjk, t) of the projective component, 
vector Hm(r-rc 0s, tk, t), where k = 1, . . . , Ndiv, 

are calculated, and tlie vector e„(tk) can be used as the direction of the noise 
10 magnetic field tlirougli utilization of the time tk in which tlie sum of tlie 
above-mentioned variances 

Z var„(.,) (58) 

or 

J Z (var„(/.)y. (59) 

\ a ^x^ y, z 

15 becomes minimum or smaller than a predetermined value (S145b). In this case, 
there is tlie possibility that the correlation functions Rx(tk), Ry(tk) and Rz(1k) at a 
certain time tk have lost tlieir original signs. Hence, it is necessary to evaluate the 
fluctuation of the projective component at each time ^ for four combinations 
[Rx(tk), Ry(tkX Rz(tk)], [Rx(tkX Ry(tk), -Rz(tk)], \RM, -Ry(tk), Rz(tk)] and [R,(tk), 
20 -Ry(tkX -Rz(tk)]. 

Further, tlie period during wliicli the signal magnetic field is OFF can 
easily be set (SI 45c, S145d) based on the sequence s(t-lk). By applying the same 
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metiiod as the first one to tlie magnetic field vector Hn(r-rc, 9c, t) measured in this 
period, the direction en(r ) of the noise magnetic field vector Hn(r, t) can be 
obtained (SI 46). In the above, use is made of the time tk of die equally spaced 
Ndiv sections of a tlie time lengtli Tdiv divided fi-om the period Tperbd, but it is also 

5 possible to use tlie time when the correlation fiinction given by Eq. (45) or (46) 
becomes maximum or when the correlation function becomes maximum and 
exceeds a predetermined value. 

(Mediod for detemiining die direction of the noise magnetic field when the 
number of noise magnetic fields is virtuallv two) 

10 Even in a case of two noise magnetic fields, when a first one of the two 

noise magnetic fields has far higlier an intensity than the second noise magnetic 
field at a first fi^equency and the second noise magnetic field has far higher an 
intensity tlian die first noise magnetic field at the second fi-equency, the directions 
of the first and second noise magnetic field can easily be calculated througli 

15 utilization of these firequency components in the measured magnetic field. 

If the first or second fi^equency is close to the fi-equency of the signal 
magnetic field, tlie directions of the noise magnetic fields can.be determined from 
the measured magnetic field having passed through a band pass filter that permits 
the passage therethrough of only frequencies near those of the signal magnetic field, 

20 by the same metliod as the fourth or fifth method for use in the case of virtually one 
noise magnetic field 

When either of the first and second frequencies does not equal to the 
frequency of die signal magnetic field, the directions of the respective noise 
magnetic fields need only to be calculated by the same method as the first one for 

25 use in the case of virtually one noise magnetic field. 
(Other embodiments) 

The present invention is also effective when the signal magnetic field 
generated by the magnetic field source is virtually axially symmetric, and the 
invention permits position determination with a smaller number of magnetic 
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sensors or by magnetic field sensing at a smaller number of positions than in die 
case of a magnetic field of low symmetiy 

Further, when only one noise magnetic field affects the measurement of tlie 
digging position and tlie tilt angle of the magnetic field source is known which is 
5 an inclination of the axis of symmetry corresponding to the axial direction of the 
signal magnetic field set in the magnetic field source with respect to tlie vertical 
direction, the projective component of the magnetic field, measured at each of two 
or more dilferent positions, on a plane perpendicular to tlie direction of die noise 
magnetic field sensed at each of the magnetic field sensing positions is calculated; 

10 the position of tlie magnetic field source and its azimuth angle that is the direction 
of the axis of symmetry in a horizontal plane can be calculated firom the above-said 
projective component. 

Further, according to this invention metiiod, when virtually one noise 
magnetic field affects the measurement of the digging position, the projective 

15 component of tlie magnetic field, measured at each of tliree or more different 
positions, on a plane perpendicular to the direction of the noise magnetic field 
sensed at each of tlie magnetic field sensing positions is calculated; tlie position of 
the magnetic field source, its tilt angle that is an inclination of tlie axis of symmetry 
corresponding to the axial direction of tlie signal magnetic field set in the magnetic 

20 field source with respect to die vertical direction, and the azimuth angle of die 
magnetic field source diat is he direction of die axis of symmetry in the horizontal 
plane can be calculated from die above-said projective component. 

Furtlier, when virtually two noise magnetic fields alone affect the 
measurement of die digging position and the tilt angle of die magnetic field source 

25 is known which is an inclination of die axis of symmetry corresponding to the axial 
direction of die signal Jiiagnetic field set in the magnetic field source widi respect 
to die vertical direction, the projective component of die magnetic field, measured 
at each of four or more different positions, on a straight line perpendicular to both 
of die direction of a first one of die two noise magnetic fields sensed at each 
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magnetic field sensing position and tlie direction of the remaining second noise 
magnetic field sensed at tlie same position is calculated; tlie position of the 
magnetic field source and its azimuth angle tliat is the direction of tlie axis of 
symmetry in the horizontal plane can be calculated from the above-said projective 
5 component. 

Further, when virtually two noise magnetic fields alone affect the 
measurement of the digging position, the projective component of the magnetic 
field, measured at each of five or more different positions, on a straight line 
perpendicular to botli of the direction of a fust one of the two noise magnetic fields 

10 sensed at each magnetic field sensing position and the direction of the remaining 
second noise magnetic field sensed at the same position is calculated; tlie position 
of the magnetic field source, its tilt angle tliat is an inclination of die axis of 
symmetry corresponding to tlie axial direction of the signal magnetic field set in 
the magnetic field source with respect to the vertical direction, and the azimuth 

15 angle that is the direction of the axis of symmetry in the horizontal plane can be 
calculated fi-oni the above-said projective component. 

Moreover, when virtually two noise magnetic fields alone affect the 
measurement of the digging position, tlie frequency component of a first one of tlie 
two noise magnetic fields, in tlie vicinity of which the remaining second noise 

20 magnetic field and tlie signal magnetic field have substantially no frequency 
components, is measured to thereby permit detection of the direction of the first 
noise magnetic field in temis of vector; and the frequency component of the 
second noise magnetic fields, in tlie vicinity of which the first noise magnetic field 
and the signal magnetic field have substantially no frequency components, is 

25 measured to thereby permit detection of tlie direction of the second noise magnetic 
field in terms of vector 

In the present invention, it is eifective to use, as a magnetic sensor, a 
three-axis magnetic sensor that senses three magnetic fields orthogonal to one 
aiiotlier at substantially the same position. 
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Tlie magnetic sensor for use in the present invention may be any kinds of 
sensors as long as they are capable to three magnetic fields orthogonal to one 
another at substantially tlie same position, but the tliree-axis magnetic sensor is 
suitable wliich senses three magnetic fields orthogonal to one another at 

5 substantially tlie same position. Altematively, it is possible that one magnetic 
sensor capable of sensing a magnetic field in only one direction is turned at the 
same position tovvaid three orthogonal directions one after anotlier to sense the 
tliree orthogonal magnetic fields. 

In carrying out the present invention, it is possible to employ such a frame 

10 12 as depicted in Fig. 32 which has magnetic sensor fixing means 11 mounted 
thereon to fix tliree-axis magnetic sensors and a tilt angle gauge 13 for detecting 
the inclination of the frame witli respect to a vertical direction. The position of 
each magnetic sensor fixing means on the frame is knovvn, and the magnetic sensor 
fixing means possesses a function of fixing tlie magnetic sensor in a predetemiined 

15 orientation to tlie frame. The magnetic sensor fixing means 11 is provided with, 
for example, tliree faces orthogonal to one another, and has a mechanism that fixes 
the magnetic sensor at a predetermined angle when a predetermined face of the 
sensor case is pressed against any one of Uie tliree faces of tlie sensor fixing means. 
One magnetic sensor or one tliree-axis magnetic sensor is fixed to tliese magnetic 

20 sensor fixing means one after anotlier to sense the magnetic fields. 

In another alternative, a plurality of magnetic sensors may be fixed in 
predetemiined orientations to tlie frame 12 at a plurality of positions to 
simultaneously sense magnetic fields at the plurality of positions. 

As described above, tlie present invention employs a fi^ime provided with a 

25 plurality of magnetic sensor fixing means each capable of removably or fixedly 
mounting a tliree-axis magnetic sensor and a tilt angle sensor capable of detecting 
the tilt angle of an orthogonal coordinate system to tlie vertical direction, tlie 
magnetic fixing means being mounted on tlie frame so that their positions and 
orientations are known; die magnetic sensor is removably or fixedly mounted on a 
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required number of magnetic sensor llxing means to sense magnetic Fields, and the 
tilt angle of tlie frame during magnetic field sensing and tlie orientation of tlie 
magnetic sensor at each magnetic sensor mounting position with respect to the 
frame are used to calculate from the magnetic field sensed at each magnetic sensor 
mounting position the sensed magnetic field, a noise magnetic field and a signal 
magnetic field as vectors in a coordinate system fixed to the ground. 

The magnetic field generating means for the signal magnetic field in the 
present invention may be a coil. The magnetic field generating means may be 
one electric wire as well, or may also be one electric wire that is straight only in the 
vicinity of the place of position determination. 

As described above, even if a buried power line, railroad tracks, or similar 
noise magnetic sources are present near a construction site, die present iiwention 
pennits liighly reliable position determination without being affected by noise 
magnetic fields generated by such noise magnetic field sources. 

The present invention is intended for measuring tlie digging position in the 
non-open-cut method of excavation, but is applicable as well to many technical 
fields that involve position determination by sensing magnetic fields. 
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